Critical Behaviour of Non-Equilibrium Phase Transitions 
to Magnetically Ordered States 
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We describe non-equilibrium phase transitions in arrays of dynamical systems with cubic nonlin- 
earity driven by multiplicative Gaussian white noise. Depending on the sign of the spatial coupling 
we observe transitions to ferromagnetic or antiferromagnetic ordered states. We discuss the phase 
diagram, the order of the transitions, and the critical behaviour. For global coupling we show ana- 
lytically that the critical exponent of the magnetization exhibits a transition from the value 1/2 to 
a non-universal behaviour depending on the ratio of noise strength to the magnitude of the spatial 
coupling. 
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Introduction. In the last decade studying arrays of 
stochastically driven nonlinear dynamical systems the 
notion of noise induced non-equilibrium phase transition 
has been established p]-p^; for a recent monograph see 
pl| . In close analogy to equilibrium phase transition one 
has order parameters and finds continuous or discontinu- 
ous transitions associated with ergodicity breaking. The 
behaviour near the transition point is characterized by 
power laws and a critical slowing down. 

In this paper we consider arrays of spatially harmon- 
ically coupled Stratonovich models which undergo 
transitions into ordered states comparable to ferromag- 
netic (FM) or antiferromagnetic (AFM) phases depend- 
ing on the sign of the coupling constant. The AFM situa- 
tion is described first in this paper for that class of mod- 
els. We determine the phase diagram and characterize 
the critical behaviour at these transitions. For the glob- 
ally coupled system we derive an analytical result for the 
critical exponent of the order parameter, i.e. the mag- 
netization. This critical exponent exhibits a hitherto not 
described transition from a value 1/2 to a non- universal 
behaviour when increasing the ratio of noise strength and 
magnitude of the spatial coupling. 

The dynamics of the individual constituents Xi at 
the lattice sites i — 1,... ,L is governed by a sys- 
tem of stochastic ordinary differential equations in the 
Stratonovich sense 
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where Af{i) denotes the set of sites interacting with site 
i. N = #A/'(i) is equal to i — 1 in the case of global 
coupling and to 2d in the case of nearest neighbour (n.n.) 
coupling on a simple cubic lattice in d dimensions. D is 
the strength of the spatial interactions. £,i{t) is a zero 
mean spatially uncorrelated Gaussian white noise with 
autocorrelation function {^i{t) £,j{t')) = a'^Sij6{t — t'), 
where cr^ is the noise strength. 

The stationary probability density Ps{xi) fulfills the 
(reduced) stationary Fokker-Planck equation M 



= 
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where {xj l^^) = / dxjXjPs{xj \xi) is the steady state con- 
ditional average of Xj, j G A/'(i), given Xi at site i. We 
denote its spatial average by 

m, = ^ {xj\x,). (3) 

Global coupling. In the case of global coupling, fluctu- 
ations of nii disappear in the limit L — > oo. We thus 
may consider rrii as a parameter and obtain except for a 
constant factor a stationary solution of (||) 



Ps{xi, mi) = \x. 



(4) 



If this expression is normalizable, the stationary proba- 
bility density Ps{xi,mi) reads 



Ps{xi,mi) = 



l/N(mi) ps{xi,mi) for Xi £ supp. 







otherwise, 



(5) 



where N{mi) — J^^^^ dx Ps{x, mi). P, lives on a support 
on which (^) is normalizable, i.e. TV is finite. 

For both D and mi nonzero the support of Ps is such 
that Dmi/xi ^ ensuring normalizability of (^). For 
TOi = normalizability requires that the exponent of the 



algebraic factor in m" 
D > a the solution {'- 



is larger than —1, i.e. D < a. For 
) is not normalizable and we have 



Ps{xi) — S{xi). The determination of is described 
below in detail. 

Varying the control parameters of the system a and 
D, or the strength of the noise cr^, one obtains the phase 
diagram shown in Fig. ^. 

We first consider D > which favours a FM order. 
In the spatially homogeneous case mi = m and for m > 
or m < the support of Pg is [0, cx)) and (— cx),0], 
respectively. 
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FIG. 2. Order parameter m vs. cr . The solid line is the 
solution of (0) compared with simulations of (jl|) for global 
(circles) and n.n. coupling in d = 1 (triangels) and d — 3 
(squares). For a < ai^' — D (a) the order parameter m 
increases monotonously with whereas for a > D (b) a pro- 
nounced minimum appears. The parameters used are a — 
in (a) and a = 1.5 in (b), D — 0.5. 



FIG. 1. Phase diagram in the {D,a) plane for global cou- 
pling. Continuous transitions toward FM or AFM states 
occur crossing the fat solid lines. The double solid line 
indicates discontinuous transitions. The critical values 
ai"' , n = 1, . . . , 4, the disordered phase, the metastable AFM* 
phase, and the inserts are explained in the text. 

All constituents have the same (statistical or temporal) 
average 

{x) — j dxxPs{x,m)^ lim — f dtx[t) — F{m)^ 

(6) 



which equals the spatial average m (ergodicity). This 
leads to the self-consistency condition 



m — {x) — F{m) 
determining m. One easily finds 



(7) 



1 for a < £). 

Eq. (|^) has a trivial stable 



(1) 



For a < a, 
solution m — Q which looses its stability at a 
is determined by the condition F'(0) = 1. It bifurcates 



.(1) 



which 



into a pair of stable solutions m 



> and TO 



— m+ — corresponding to a continuous transition 
from a paramagnetic to a FM situation. Choosing to = 
TO+ for instance, the stationary probability distribution 
of the corresponding ergodic component is Ps{x, to+), cf. 
Eq. (|). In the FM region, for ai^'' < a < ai^^ = D the 
magnetization to — (x) increases monotonously with ct^, 

(2) 

whereas for a > ac there is a nonmonotonous behaviour, 
cf. Fig. |. 

As a function of D, the magnetization m increases 
continuously from zero when increasing D from zero for 
ai^^ < a < 0, whereas the transition is discontinuous for 
a > as shown in Fig. |^, cf. also Q. 



Within the FM region, a metastable antiferro- 
magnetic solution (AFM*) exists besides the stable FM 
solution for a > a^c^\ see Fig. ^ The critical value 
ai^^ — 3/2D + (7^/2 is obtained for weak noise from an 
extremal approximation for to = (x) in the spirit of Q . 
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FIG. 3. Order parameter m vs. D. Lines and symbols have 
the same meaning as in Fig. |2| For — = ai^' < a < 
the transition is continuous (a), whereas for < a it is dis- 
continuous (b). The parameters used are a = —0.5 in (a) and 
a = 1 in (b), = 2. 



(2) 

For D < the situation is different. For a < Uc 

(2) (3) 

we have Ps{x) — S{x). In the range ac < a < Uc — 
D + 0-'^ /2 one finds m — and the stationary probability 
density Ps{x) lives on (—00,00), we call this the disor- 
dered phase. For a > Oc^'' the stationary solution (^) 
is normalizable only for to 7^ 0, for to > or to < 
the support is (— oo,0] or [0, 00), respectively. We define 
two subsystems labelled by -I- and — for which the aver- 
ages (xi) have -I- or — sign, respectively. For global cou- 
pling AFM order implies to± =pO in the limit L — > 00. 
Therefore, the mean values (xi) — {x±) are given by 

{x±) = -{x^)^ dxxP,{x,TO)^±F{±0), (9) 

Jo 

where Pg is taken from (|^). 
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Nearest neighbour coupling. For n.n. coupling on a cubic 
lattice a mean field approximation is obtained in a sim- 
ilar way replacing the spatial average over the 2d near- 
est neighbours as rui = l/i'2d)J2jeJ\f(i)(^j\^i) ~ (2^*)- 
The FM case, D > 0, is formally the same as for global 
coupling but Eq. holds only approximately. In the 
AFM case, D < 0, one should take into account that 
now the two subsystems -I- and — correspond to differ- 
ent Neel sublattices A and B, respectively, and all the 
nearest neighbours of a given lattice site belong to the 
complementary sublattice. Self-consistency requires 



m± 



(10) 



For £> < system (0) is invariant under the transfor- 



for i e A, Xi 



for j e B, 



mation [p_5[ xi 
D —> —D, a ^ a — 2D. This implies that properties of 
the AFM phase for spatial coupling D ^ —D' < can 
be inferred from properties of the FM phase for spatial 



coupHng strength £>', cf. 

.(3) 



If 



|]. For instance, ai^^ — — 
transforms into af = 2D — a^/2. The phase diagram 
for n.n. coupling is shown in Fig. 0. 
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FIG. 4. Phase diagram in the (D, a) plane for n.n. cou- 
pling. For D > the situation is the same as for global 
coupling except that the metastable AFM* phase is absent 
here. For the case D < see text. 



Critical behaviour. Varying the control parameters a 
and one observes continuous transitions from zero to 
nonzero values of m with a characteristic power law be- 
haviour near the critical values of the control parameters. 
To analyze the critical behaviour it is useful to write the 
self-consistency equations and ( p^ ) in compact form 
as 



2\D\ fd\nl{m) 
dm 



(11) 



where 



In the limits cr —> or £) ^ cxd this integral can be 
evaluated by the Laplace method, cf. e.g. 0. Inserting 
the results in (|l^) for small m, one obtains the power laws 
m {a + + \D\ - D)^/'^ for cr -> and m {a + 
^2^2)1/2 for D ^ 00 with the critical exponent /3 — 1/2, 
cf. also (2|,|. 

For finite values of a and D the scaling behaviour of 
/(m) can be evaluated for small m with the result ||T^ 

/(to) ~ to2(-^^)/-' (1 -I- C7iTO-2(--^)/-^ + C2to2), (13) 

where e = a — Qc- The critical value Oc is — ct^/2 for the 
FM case and 2D - a^/2 for the AFM case. Inserting (|l|) 
in Eq. (^ij) we obtain for small m and in lowest order of 
e the power law 

/5 = sup{l/2,aV(2|D|)}, 



(14) 



logarithmic corrections are easily computed. Obviously, 
the critical exponents are the same varying a or ct^, i.e. 
Pa — Pa — P using notations from For models where 
the cubic nonlinear! ty in (|l|) is replaced by x^^^ , p > 0, 
in (^ the value 1/2 is replaced by 1/p 

Fig. 1^ compares the magnetization m{a) and the criti- 
cal exponent /3 obtained from the analytical results with 
simulations for both global and n.n. coupling [Q. For 
global coupling, simulations for systems of size 10'^ are 
already very close to the results for the infinite system. 
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FIG. 5. Critical behaviour. The Figure shows the order 

parameter m vs. control parameter a for different values of D 

2 



0.1, i.e. 



-0.05) and the critical exponent (3a vs. 



a /{2D). (a) and (b) refer to global couphng, (c) and (d) to 
n.n. coupling in d = 3 , respectively. In (a) and (c) the solid 
line is the numerical solution of the self-consistency equation 
(0). In (b) and (d) the solid line is the analytical result (14); 
the + symbols represent the numerical solution of Eq. (LI). 
Circles and squares result from simulations of Eq. (|l|); the 
triangle in (d) is the numerical result of Error bars are 
partially smaller than the symbol size. 
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Conclusion. For the globally coupled model we found 
analytically a transition of the critical exponent /3 from 
a value 1/2, which reflects the order of the nonlinearity 
and is independent of the strength of noise and spa- 
tial coupling D, to a non-universal behaviour, depending 
on (T^ and D independent of the order of the nonlinear- 
ity. This differs from the value /? = 1 proposed for the 
continuous version of the model in ||^ . 

If the noise is not too strong, the 'mean field' results 
describe the critical behaviour for n.n. coupling observed 
in our simulations very well, cf. Fig. |^d. Also the nu- 
merical result /3 w 1 obtained by Genovese and Mufioz 
1^ near a = — 1, ct^ = 2 for Z? = 1 (their 'weak noise 
phase') is in accord with our analytical result (|lj). For 
stronger noise, simulations for n.n. coupling may differ 
considerably from the mean field prediction. 

For models with a nonlinearity 2;^"'"^, the value 1/p of 
P is in general different from the value 1/2 characteristic 
for models with only additive noise. In this case, one 
has to expect interesting crossover phenomena for models 
with additive and multiplicative noise when changing the 
relative strength of the noises. 

This study was partially supported by the DFG with 
Grant BE 1417/3. 
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